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Theory of quantum electrodynamics in three spatial-time dimension is applied to the two- 
dimensional 5=1/2 quantum Heisenberg antiferromagnet in order to investigate a doped hole 
in high-temperature superconductors. Strong coupling analysis of the U(l) gauge field interaction 
is carried out to describe spectral broadening observed in the undoped compounds. It is found that 
the fermionic quasiparticle spectrum is of Gaussian form with the width about 3J, with J being the 
superexchange interaction energy. The energy shift of the spectrum is on the order of the quasipar- 
ticle band width, which suggests that the system is in the strong coupling regime with respect to 
the gauge field interaction describing the phase fluctuations about the staggered flux state. 
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I. INTRODUCTION 



One of the most fundamental questions about high-temperature superconductivity is how to describe doped holes 
introduced in the Cu02 plane. The simplest way to approach this problem is to investigate the single hole doped 
j_i ■ system. Experimentally, the excitation spectrum associated with a single hole is observed by angle resolved pho- 
toemission spectroscopy (ARPES) in the undoped compounds such as S^CuC^Clzii^ and Ca2Cu02Cl2. 3 Although 
(_J . the spectrum is not sharp but quite broad, whose width is ranging from O.leV to 0.5eV, the experiments show that 
the trace of the peak indicates a dispersion whose maxima are at (±7r/2, ±7r/2). The dispersion near these points is 
quadratic and almost isotropic. The band width is 2.2 J ~ 270meV with J the superexchange interaction. This band 
width is much smaller than the band-structure estimation of 8t ~ 2.8eV. Furthermore, the observed spectra are not 
described by a conventional Lorentzian form but described by a Gaussian form.— These observations suggest that the 
quasiparticle excitations in the undoped cuprates are quite different from conventional Fermi liquid quasiparticles. 
£j | In the slave- fermion theory of the t-J model, the single hole system has been analyzed by the self-consistent Born 
'— approximation^ The effect of the spin- wave excitations is included in the self-energy in a self-consistent manner 
with omitting vertex corrections. The resulting hole dispersion has minima at (±7r/2, ±7r/2), and the band width is 
scaled by J. The dispersion along (0, 0) to (tt, 7t) is in good agreement with the experiments. Quantum Monte Carlo 
simulations based on a model, in which canonically transformed spinless fermions propagate with antiferromagnetic 
\ spin correlation background, like the slave-fermion formalism, are consistent with this result^ However, the dispersion 
I/"") ■ along (tt, 0) to (0, tt) is much smaller than that in the experiments. This discrepancy is improved by including the 
next nearest neighbor and the third nearest neighbor hopping terms' This suggests that within this approach the 
quadratic behavior near the (7r/2, 7r/2) point along (0, tt) to (tt, 0) has a different origin from that along (0, 0) to (tt, tt). 
Furthermore, it turns out that the damping effect due to the coupling with the spin-wave modes does not lead to 
a broad line shaped Recently Mischenko and Nagaosa studied a coupling to an optical phonon^ They numerically 
summed over Feynman diagrams including vertex corrections for phonons. It was argued that the coupling is in the 
, strong coupling regime, and so the quasiparticle spectrum becomes broad. In this scenario the most enigmatic feature 
of the hole spectral broadening in the undoped compounds is associated with phonon effects. The dominant role is not 
r> ] played by the antiferromagnetic correlations which is believed to be essential for the mechanism of high-temperature 
jj^ ■ superconductivity. 

Here I take a different approach. I consider the staggered flux state proposed in the literatur e) 10 ' 1 from a mean 
field theory of the 5=1/2 antiferromagnetic Heisenberg model based on a fermionic representation of the spins. The 
dispersion of the quasiparticle in the staggered flux phase is in good agreement with the experimentally obtained 
dispersion as pointed out by Laughlin. 12 Including phase fluctuations about the mean field, the effective theory is 
described by quantum electrodynamics in three spatial-time dimension, which is called QED3. At mean field level, 
the fermions are massless. By including the effect of the gauge field, the mass of the Dirac fermions is induced^ This 
mass is associated with the staggered magnetization. The presence of the mass term is also suggested by a variational 
Monte Carlo approach.— The quadratic dispersion around (tt/2,tt/2) observed in the experiments is consistent with 
the massive quasiparticle spectrum. Furthermore, the quasiparticle dispersion is isotropic at (tt/2,tt/2). 

The purpose of this paper is to argue that the coupling of the fermions with the phase fluctuations leads to a broad 
Gaussian spectrum. The QED3 action with the mass term is analyzed by performing a canonical transformation. 
The spectral function is obtained by calculating the Green's function of the Dirac fermion which is associated with a 
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single quasiparticle propagation. It is shown that the spectral function shows a broad Gaussian peak whose width is 
on the order of J. 

The rest of the paper is organized as follows. In SecHH we describe the QED3 formalism of the staggered flux 
state. After taking the transverse gauge, a canonical transformation is applied. The quasiparticle Green's function 
is computed in Sec lIIIl The spectral function is obtained with including vertex functions arising from random phase 
approximation about the instantaneous longitudinal gauge field interaction. Implications of the result is discussed in 

EY1 

II. QED 3 THEORY OF THE STAGGERED FLUX STATE 

For the description of the S = 1/2 two-dimensional quantum Heisenberg antiferromagnet, I take the following 
QED3 action in the real time formalism as the effective theory: 



S= d 



-;3 ■ 



i\> (x) [i^ (d„ - ia^) - ma 3 ] ip (x) - -^UuF" 



(2.1) 



where the Dirac fermion fields, tp(x), consist of four component associated with even and odd sites and two independent 
nodes. Due to the spin degrees of freedom there are two species of tp(x). (For the derivation, see Appendix [A"l ) 
Hereafter the spin index for ip(x) is suppressed because the spin degrees of freedom does not play an important role 
in the following calculation. The action (|2.ip describes low- lying excitations around (±7r/2, ±7r/2) because the above 
continuum model was derived by taking the continuum limit at those points. Note that the theory is particle- hole 
symmetric. Therefore, the quasiparticle properties are identical to the quasihole properties. 

For the gauge, I take the transverse gauge: V • a = 0. In this gauge, the interaction between the fermions mediated 
by the longitudinal part of the gauge field is instantaneous as in the conventional electromagnetic field formulation: 



S = I d 3 xt/j (x) [ij°dt + i"f (dj + iaj) — 7710-3] ^ {%) 

^Jd 3 xJ d 2 r' [p (r, t) - po] [p (r', t) - p ] In |r - r'| 



4tt 
1 

'2^2 



3„ 



(5 t a) 2 - (V x a) 2 , (2.2) 



where the background gauge charge, —e a po, comes from the constraint on the fermion number to represent the 
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spin 1/2. In three spatial-time dimension, the "Coulomb" interaction is described by V (r) = — |Mnr. Under the 
transverse gauge, the vector potential is represented by 

a x (q) = a(q), a y (q) = —a(q), (2.3) 

q q 

in the momentum space. Quantizing the transverse gauge field, the Hamiltonian is 



H 



J d 2 riP (r) (-<y A,- + ma 3 ) 1> (r) + J d 2 r ^ (r) J (q xJv - lx q y ) ^ (b q + bl q ) V> (r) 

+ l - J d 2 r J d 2 r'V{\v~r'\)[p{v)-p ] [p(r')-p ] 
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In order to investigate the strong coupling effects, I perform the following canonical transformation^: 

H = e s He- s , (2.5) 

where 

s = Jd 2 r (b q bl q ) (r) M q (r) ^ (r). (2.6) 
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(The area of the system is set to unity.) This is a unitary transformation if M^_ q (r) = M q (r). The function M q (r) is 
chosen so that the interaction term between the Dirac fermions and the gauge field is cancelled by [s,H]: 



M q (r) = -JJ^e***^ {q xlv - lx q y ) . (2.7) 

Under this canonical transformation, the fermion fields transform as 

e s il>{v)e~ s = e x(r) V(r). (2.8) 

The function X (r) is, 

III. QUASIPARTICLE GREEN'S FUNCTION 

Now I compute the time-ordered Green's function: 

G(r,t) = -i(Tip(r, t)^ (0, 0)) 

= -i(Te x ^(r > t)^(0,0)e- x ^)_. (3.1) 

This Green's function has a matrix form of 4 x 4. But the matrix is divided into two blocks in which each part 
describes either the Dirac fermion fields near (tt/2, tt/2) or the Dirac fermion fields near (— tt/2, tt/2). It is enough to 
focus on one of them because two components are decoupled as far as long-wave length gauge field fluctuations are 
concerned. The superscript (1) is used to denote the former. Diagonalizing the factor with 7 matrices, X(r) is 



X« (r) = i£ (b q e^ + b\e-^ r ) ^U q r 3 Ul (3.2) 

where 

U q =±( ( J. w -(&-^)/9V (3.3) 

The first term in Eq. (|2.4p . which does not change its form by the canonical transformation, is diagonalized by a 
unitary transformation as well. Finally, by making use of the following formula, 

e Ab^+Bb e Cb i + Db^ _ e ±(AB+CD+2BC) e (A+C)(B+D)n(uj) ^ ^ 

for bosons, Eq. (|3.ip is 

G (1) (r, t)=-iJ2 e lW r e- lEkt R k exp [-K(r, *)] , (3.5) 

k 

at T = 0. The retarded Green's function has the same form for t > 0. Here, 



* 11 [4 V gVA: 2 + m 2 



(3.6) 



X ( r '*) = E^3 (l-e-^ 4 e^ r ). (3.7) 



Performing the Fourier transform, the spectral function A^ (k, uj) is given by A^(k,uj) = — ^ImG^ (k, u). 

So far the bare vertex function has been used for the computation. However, in the long-wave length limit taking 
the bare vertex is not appropriate as manifestly seen by infrared divergence in K(v, t). For the vertex part, random 



4 



phase approximation is applied with respect to the longitudinal interaction term. The bare vertex is reduced by the 
factor of 1/(1 — TT q ), where 

n q = v q V Tr [G kl0 G k+ql0 } =-^\ + 0(1), (3.8) 
k 11 q 

with v q = e 2 a /q 2 . For this computation, it is convenient to use the Euclidean formalism because the main contribution 
comes from ir q with q = (q, 0), where the Minkowskiian formalism leads to the same result. 
Including the vertex correction the function K(r, t) is, 



^*) = ff4(r^) 2[1 - e ^ tJo(9r)] 



(3.9) 



with Jq{x) the zero-th order of the Bessel function of the first kind. The ultraviolet cutoff A introduced here because 
the wave length of fluctuations is larger than the lattice constant. The integrand is expanded with respect to q, before 
the integration. The result is 

K(r,t) ~ iE s t+ -A 2 r 2 + -A 2 t 2 , (3.10) 
8 4 



where 
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A 2 = ^. (3.12) 
2tt 

The first term in Eq. (l3.9|) represents the energy shift, E s . The shape is changed to a broad Gaussian form by the 
subsequent terms as shown below. For the factor R k , the analytic expression was not obtained. From a numerical 
computation I found that R k is linear in k at < k < k c , with fc c ~ 1, and reaches a saturated value of 0.2 for 
k > k . To approximate R k , I took an approximate form of R k ~ 0.2fc. In computing the Fourier transform, it is 
useful to note that the integration with respect to r and that with respect to t are carried out separately. For t«A, 
I obtained 



A(X)n ^ °- 10 
A (1, (k, (J) ~ exp 



(u-E k - E s 



A 2 



(3.13) 



Therefore, the Dirac fermion energy spectrum is shifted by E s and is a Gaussian form with the width of A. For 
the estimation of these values, the spin wave velocity is assumed to be, c sw = 1.6 J. The mass term is evaluated as 
m ~ 1.3 J from the ARPES experiments by fitting the dispersion near (tt/2,tt/2). The gauge charge is simply taken 
from the factor of the Maxwellian term obtained by integrating out the Dirac fermions: e 2 — inm. (If the same 
calculation is carried out for Dirac fermions with k > ko, then the gauge charge value is e 2 ~ 37r 2 fco/4, for fco m. 
Therefore, the above choice is the minimum value for the gauge coupling.) Recalling the fact that c sw is taken unity 
in the above calculation, I found E s ~ 2 J and A ~ 3 J. This value of A is consistent with the above assumption 
about k. The weight of the Gaussian spectrum is ~ 0.05 for this value of A. Because of the spin degrees of freedom 
and the degenerate nodes, the weight is ~ 0.2 in total. From the value of E s , one can get insight about the strength 
of the coupling. If the system is in the (weak) strong coupling, the value of E s is expected to be large (small). The 
fact that E s is on the order of the band width suggests that the coupling is in the strong coupling regime. 



IV. DISCUSSION 



In this paper, the spectral function of the quasiparticle in the staggered flux state with phase fluctuations was 
computed within the effective theory described by the QED3. It was shown that the quasiparticle spectra become 
a broad Gaussian form with an energy shift due to the gauge field interaction. The estimated spectrum width is 
consistent with the experiments. The analysis suggests that the coupling to the gauge field is in the strong coupling 
regime. 
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Since the model is based on the continuum approximation, the result is applicable to the quasiparticle excitations 
near (±7r/2, ±7r/2). However, the result is extended to other k points by formulating the theory on the lattice. Such 
a model is useful to study the change of the width of the spectra away from (±7r/2, ±7r/2). 

As for the vertex correction, random phase approximation is applied with respect to the interaction arising from 
the longitudinal component of the gauge field. Of course, this is not the full vertex correction. In the long- wave length 
limit, there are other intermediate processes. However, it is expected that dominant contribution is covered by the 
above vertex correction because the longitudinal component plays a major role in screening the gauge charge. 

Finally, let me comment on vanishing quasiparticle spectra observed in the experiments^ along the line from 
(ir/2, 7r/2) to (it, 7r) and that from (it, 0) to (ir, it). A similar behavior is also observed in the pseudogap phased that 
is, only a part of the Fermi surface is observed as an arc shaped One might expect that damping effect coming from 
the coupling to the gauge field leads to the suppression of the quasiparticle peaks. However, it turns out that the 
inclusion of a slight hopping term in the staggered flux state leads to the vanishing of the spectrum in the second 
magnetic Brillouin zone. I will discuss this matter in a future publication. 
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APPENDIX A: DERIVATION OF QED 3 ACTION 

In this appendix, I derive the QED3 action as the effective theory for the S = 1/2 two-dimensional Heisenberg 
antiferromagnet : 

H = jJ^S.-Sj. (Al) 

A fermion representation is introduced for the spin 1/2: Sj^ = fj a <J a j3fjp/2 (/j = x,y,z). are the Pauli spin 

matrices. These fermions need to satisfy the constraint, ^2 a fj a fja = 1- Introducing Lagrange multipliers to take 
into accout the constraint, the Hamiltonian is, 

h = -\ j£ fifiafy* + e ^ {flf* - 2s ) > ( A2 ) 

(i-'j) 3 

up to a constant term. The mean field taken in the 7r-flux state is Xij — ^/ja/iaVh, by choosing a suitable gauge.— 

Since the system is homogeneous, uniform Xij an d Aj are assumed: xi = Xj+£,j, X2 — Xj,j-yi X3 = Xj-x,j, an d 
Xa = Xj,j+yi with j residing at an even site. Numerically solving the mean field equations for Xj (j = 1> 2, 3, 4), with 
setting A = 0, the 7r-flux stated is found in which X1X2X3X4/IX1X2X3X4I = — 1- 
The quasiparticle energy dispersion in the staggered flux state is 

E k = ±^\ X ie- lk * + X *2^ + X se ikx + xl^H (A3) 

\Ek\ has minima at (±7r/2, ±7r/2), and around these points the energy dispersion has the relativistic form. Introducing 
the even and odd site fields, f e ka = (fka + fk+Q,a)/V2 and f Q ka = (fka~ fk+Q,a)/V2 with Q = (tt,tt), and expanding 
Eh around (±7r/2, 7r/2), the Hamiltonian is rewritten as, 

TT ^ ,VV ft ft \( -X\k x +X2ky \( felka \ 

k ' \ -Xlkx + X* 2 ky J { fol ka J 

,/VVft ft \f X*kx+X2ky\ ( fe2ka\ , M] 

+ J 2^ V h2k. a lo2 ka >\xik x + X2\ky \ f o2ka J ■ [A ^> 

k x 
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Here the indices 1 and 2 are introduced to denote the fields around (7r/2,7r/2) and those around (— tt/2, 7r/2). The 
summation with respect to k is taken over the magnetic Brillouin zone: \k x ± k y \ < tt. Choosing \i — X3 — i\x\ an d 

X2 = Xi = \x\, and setting V>L = ^f\ lka , f olka , f o2ka , / e f 2ta ) , the action is, in the continuum limit, 

S = { d 3 xlp (x) i^d^ (x) , (A5) 
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where ip — ip'j and the 7 matrices are 

7° = 

Phase fluctuations about the staggered flux mean field state are included by the U(l) gauge field, a„, by replacing 
dfj, with — ia^. Integrating out the high-energy Dirac fermion fields with k > kg leads to the dynamics of the 
gauge field, which has the Maxwellian form. Numerically solving the Schwinger-Dyson equation, it is found that the 
self-energy has a non-zero massJ^ Physically this mass is associated with the presence of the staggered magnetization 
which is absent at the mean field level. From the variational Monte Carlo approach with the finite mass m, it is 
shown that the mean field energy improves by the inclusion of mi 19 ' 20 The action (|2.1|) is obtained by including the 
mass term arising from the staggered magnetization. 



* morinari@yukawa.kyoto-u.ac.jp 

1 B. O. Wells, Z. X. Shen, A. Matsuura, D. M. King, M. A. Kastner, M. Greven, and R. J. Birgeneau, Phys. Rev. Lett. 74, 
964 (1995). 

2 S. LaRosa, I. Vobornik, F. Zwick, H. Berger, M. Grioni, G. Margaritondo, R. J. Kelley, M. Onellion, and A. Chubukov, 
Phys. Rev. B 56, R525 (1997). 

3 F. Ronning, C. Kim, D. L. Feng, D. S. Marshall, A. G. Loeser, L. L. Miller, J. N. Eckstein, I. Bozovic, and Z. X. Shen, 
Science 282, 2067 (1998). 

4 K. M. Shen, F. Ronning, D. H. Lu, W. S. Lee, N. J. C. Ingle, W. Meevasana, F. Baumberger, A. Damascelli, N. P. Armitage, 
L. L. Miller, et al., Phys. Rev. Lett. 93, 267002 (2004). 

5 C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 6880 (1989). 

6 M. Brunner, F. F. Assaad, and A. Muramatsu, Phys. Rev. B 62, 15480 (2000). 

7 T. Tohyama and S. Maekawa, Supercond. Sci. Technol. 13, R17 (2000). 

8 J. Bala, A. M. Oles, and J. Zaanen, Phys. Rev. B 52, 4597 (1995). 

9 A. S. Mishchenko and N. Nagaosa, Phys. Rev. Lett. 93, 036402 (2004). 

I. Affleck and J. B. Marston, Phys. Rev. B 37, 3774 (1988). 

1 X. G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996). 

2 R. B. Laughlin, Phys. Rev. Lett. 79, 1726 (1997). 

3 D. H. Kim and P. A. Lee, Ann. Phys. (N. Y.) 272, 130 (1999). 

4 T. C. Hsu, Phys. Rev. B 41, 11379 (1990). 

G. Mahan, Many-Particle Physics (Plenum publishers, New York, 2000), 3rd ed. 

6 A. Damascelli, Z. Hussain, and Z. X. Shen, Rev. Mod. Phys. 75, 473 (2003). 

7 M. R. Norman, H. Ding, M. Randeria, J. C. Campuzano, T. Yokoya, T. Takeuchi, T. Takahashi, T. Mochiku, K. Kadowaki, 
P. Guptasarma, et al., Nature 392, 157 (1998). 

8 P. A. Lee, N. Nagaosa, and X. G. Wen, Rev. Mod. Phys. 78, 17 (2006). 

9 C. Gros, Ann. Phys. (N. Y.) 189, 53 (1989). 

!0 T. K. Lee and S. P. Feng, Phys. Rev. B 38, 11809 (1988). 



